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Abstract 

We review a simple model of closed economy, where the economic agents make 
money transactions and a saving criterion is present. We observe the Gibbs distri- 
bution for zero saving propensity, and non-Gibbs distributions otherwise. While the 
exact solution in the case of zero saving propensity is already known to be given 
by the Gibbs distribution, here we provide the explicit analytical form of the equi- 
librium distribution for the general case of nonzero saving propensity. We verify it 
through comparison with numerical data and show that it can be cast in the form 
of a Poisson distribution. 
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1 Introduction 



It is known that the higher end of the distribution of income f(m) follows 
the Pareto law [1], f(m) oc m _1_a , where m is the income (money) and the 
exponent a has a value in the interval 1 and 2 [2,3,4,5]. An explanation of 
the Pareto law, in terms of the laws regulating the system micro-dynamics, 
should take into account its basic constituents, i. e. the trading agents, as well 
as the criteria used to carry out the economic transactions. Several studies have 
been made to provide an explanation (see Ref. [6] for a brief summary and 
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more references). In this respect, it is of general interest to study some simple 
systems of closed economy, which can be either solved exactly or simulated 
numerically, in order to investigate the relation between the micro-dynamics 
and the resulting macroscopic money distribution [7,8,9,10,11,12]. In this pa- 
per we consider the generalization of a simple model of money conserving 
economy, realized by introducing a criterion of saving in the transaction law, 
through the saving propensity A. We study numerically its asymptotic money 
distribution as a function of the model parameters. We show that it is not a 
Gibbs distribution and, by direct comparison with numerical data, that the 
corresponding analytical solution has the form of a Poisson distribution. 



2 Model 



In the simple model considered [8], N agents can exchange money in pairs 
between themselves. For the sake of simplicity we assume that all the agents 
are initially assigned the same money amount mo, despite this condition is 
not restrictive for the following results. Agents are then let to interact. At 
every "time step", a pair is randomly chosen and the transaction takes 
place. During the transaction, the agent money amounts rrii and rrij undergo a 
variation, rrii — > w! i and rrij — > m'j. Money is assumed to be conserved during 
the transaction, so that 

rrii + TUj = m' i + m'j . (1) 



In this basic model, m- and m'j are obtained through a random reassignment 
of the total money (rrii + m j), 



m i = e(m i + rrij) , 

m' = (1 - e)(rrii + rrij) , 



(2) 



where e is a random number, extracted from a uniform distribution in the 
interval (0,1). Notice that this model of dynamics, as well as its variations 
considered in the following, ensures that agents have no debts after the trans- 
action, i. e. they are always left with a money amount m > 0. It can be shown 
that, merely as a consequence of the conservation law (1), the system relaxes 
toward an equilibrium state characterized by a Gibbs distribution [7,8,9], 

f(m) = f3exp (-(3m) , (3) 



where (3 — 1/ (m) represents the inverse average money and (m) = J2i m i/N = 
rriQ. This means that after relaxation, the majority of the agents has a very 
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small amount of money, while the number of richest agents - e. g. those with 
m larger than a given value m', as well as the fraction of the total money 
they own, exponentially decreases with m' . The Gibbs distribution (3) has 
been shown to represent a robust equilibrium state, reached independently 
of the initial conditions also in generalized models, such as those involving 
multi-agent transactions. 

However, if a saving criterion is introduced [7,9], i.e. agents save a fraction 
A - the saving propensity - of the money they have before the transaction 
is made, the shape of the equilibrium distribution changes dramatically. The 
conservation equation (1) still holds, but the money to be shared in a trans- 
action between the i-th and the j-th agent is now (1 — A)(mj + rrij). Then 
Eqs. (2) are thus modified, 



m\ = Xrrii + e(l — A) (to* + rrij) , 

m'j — Xrrij + (1 - e)(l - A)(m; + rrij) . (4) 
These equations can also be rewritten in the following way, 



m' i = m i + Am , 
m'j = mj — Am , 

Am = (1 - \)[errij - (1 - e)m,] , (5) 

which clearly shows how money is conserved during the transaction. 

We performed numerical simulations, for various values of A, of a system with 
N = 500 agents. In each simulation a sufficient number of transactions, as far 
as 10 7 , depending on the value of A, was used in order to reach equilibrium. 
The final equilibrium distributions for a given A, obtained by averaging over 
1000 different runs, are shown in Fig. 1. 



3 Fitting 



The exact solution for the case A = is known to be given by the Gibbs 
distribution, Eq. (3). Here we give the corresponding exact solution for an 
generic value of A, with < A < 1. This solution was found by fitting the 
results of numerical simulations and it turns out to fit extremely well all data. 

It is convenient to introduce the reduced variable 
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Fig. 1. Equilibrium money distributions for different values of the saving propensity 
A, in the closed economy model defined by Eqs. (5). The continuous curves are the 
fitting functions, defined in Eq. (8). 
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Fig. 2. As in Fig. (1), but on a double logarithmic scale, 
the agent money in units of the average money (to), and the parameter 

3A 

n( A) = l + — . 



(7) 



We found that the money distributions, for arbitrary values of A, are well 
fitted by the function 

P(x) = a n x n ~ l exp (—nx) , (8) 
where x and n are defined in Eqs. (6) and (7), respectively 1 . Using normal- 



1 An excellent fitting is also obtained if the variable in the exponential is raised to a 
power c, exp(— nx) — > exp(— nx c ), where c is an additional parameter and the value 
of c is close to one for all values of A. Here we assume c = 1, since the corresponding 
fitting is good. 
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ization conditions, the prefactor is easily shown to be given by 



n" 



a n = 



r(n) ' 



(9) 



where T(n) is the Gamma function. 

The fitting curves for the distribution (continuous lines) are compared with 
the numerical data in Fig. 1. The fitting describes the distribution also at 
large values of x, as shown by the logarithmic plots in Fig. 2. The numerical 
values of the parameters a n and n are compared with the respective fitting 
functions (9) and (7) in Fig. 3. The distribution function (8) still contains an 
exponential factor exp(— nx), similar to that of the Gibbs distribution, but the 
average value is now rescaled by n. The power x n_1 qualitatively changes the 
Gibbs distribution into a curve with a maximum at x > 0, i.e. with a mode 
different from zero. 




Fig. 3. The parameters a n (left) and n (right) versus A obtained from numerical 
data (dots) and the corresponding analytical formulas (continuous curves) given by 
Eqs. (9) and (7), respectively. 



It is to be noticed that by introducing the rescaled variable 

7-7-7- , (10) 

{m) jn 



and the corresponding probability density P n (x n ) = dF(x)/dx n = P(x)/n, 
where F(x) is the cumulative function, and using the explicit expression of 
a n , Eq. (9), the distribution (8) becomes 

P„(x n ) = — 1 <~ 1 exp(-o; n ) , (II) 
r(n - 1) 



which reduces to the Poisson distribution for integer values of n. 
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4 Conclusions and discussion 



We have studied a generalization of the simple closed economy model, in 
which a random reassignment of the agent money takes place, by introducing 
a saving propensity A > 0. We have empirically obtained the corresponding 
exact analytical solution from a fitting of the numerical data. The distribution 
naturally lends itself to be interpreted as a Poisson distribution P(n, x n ) for 
the reduced variable x n = m/((m) /n). The parameter n= 1 + 3A/(1 — A) is in 
principle continuous but it can vary between 1 and oo when A varies between 
and 1. This result raises the problem of a more rigorous derivation of the 
solution as well as of a deeper physical interpretation of the result. 
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